Thermalization process of a photo-generated plasma in semiconductors 
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The kinetics of ultra-fast processes which leads to the thermalization condition of a photo-excited 
plasma in semiconductor systems is studied theoretically. We analyze the time evolution of a carrier 
population generated by a finite optical pulse, from the beginning of the pulse until the time in 
which the carrier population reaches a quasi-equilibrium condition. We calculate the energy fluxes 
caused by the main interaction mechanisms along the different stages the system passes through. 
Our analysis is done by using a set of non-linear rate equations which govern the time evolution of 
the carrier population in the energy space. We consider the main interaction mechanisms, including 
dynamic screening and phonon population effects. 

Keywords: photo-excited plasma; thermalization; ultrafast processes in semiconductors 

Se estudia la cinetica de los procesos ultra rapidos que llevan a la condicion de termalizacion 
de un plasma fotoexcitado en sistemas semiconductores. Analizamos la evolution temporal de 
una poblacion generada por un pulso optico finito, desde el comienzo del pulso hasta el tiempo 
en el que la poblacion alcanza una condicion de cuasi equilibrio. Calculamos los flujos de energfa 
causados por los mecanismos principales de interaction a lo largo de las diferentes etapas por las que 
pasa el sistema. Hacemos nuestro analisis usando un conjunto no lineal de ecuaciones de razon de 
cambio que gobiernan la evolution temporal de la poblacion de portadores en el espacio de energfas. 
Consideramos los mecanismos principales de interaction, incluyendo el apantallamiento dinamico y 
efectos de poblacion de fonones. 

Descriptores: plasma fotoexcitado; termalizacion; procesos ultra rapidos en semiconductores 
PACS numbers: 72.20.Jv, 72.20.Dp, 78.20.Bh, 78.47.+p 
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I. INTRODUCTION 

The relaxation processes of photo-generated plasma 
systems in semiconductors exhibit two characteristic 
stages. The first of them, commonly called the ther- 
malization process, is mainly governed by the rapid in- 
teractions, namely, the electron-electron (e-e) and the 
electron-optical phonon interactions. In this stage the 
carrier distribution function (CDF) is far from equilib- 
rium and theoretical approximations for this, based on 
small displacements from the equilibrium are poor ap- 
proaches to describe the kinetics of these processes. This 
stage ends when the interaction mechanisms random- 
ize the energy and momentum in the carrier population. 
This allows the CDF to reach a condition in which it is 
possible to define an effective temperature for the carrier 
system, i.e., the CDF acquires a shape similar to that of 
an equilibrium one. The second stage of the relaxation is 
the so called, cooling process, and is mainly ruled by the 
slow interactions in the system, namely, electron-phonon 
(e-ph) scattering and recombination. This process has 
been extensively studied since the pioneering works of the 
late 60's and 70's decadesQ, 0. In this stage, the time 
evolution of the quasi-equilibrium CDF can be described 
by means of simplified evolution equations for time de- 
pendent effective temperature and chemical potential 
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There have been in the literature reported theoretical 
and experimental studies on the thermalization process, 

however it is still not well understood @,HH HQ. ThiS 
is because the thermalization process is a stage domi- 
nated by transient effects in a far from equilibrium sys- 
tem. The understanding of the transient processes oc- 
curring in photo-generated carrier populations is of great 
relevance because it would allow a deeper physical insight 
on the dynamical effects of interaction mechanisms upon 
observable properties of the system. In addition, it could 
provide information, as well as a theoretical framework, 
to investigate some particular ultra-fast phenomena like 
the kinetics of thermo-transport and the kinetics of the 
coherent control of quantum states in mesoscopic semi- 
conductor systems fcjlfjlfioj. 

Thermalization depends mainly on the efficiency of the 
e-e scattering to redistribute the excess energy given to 
the system by the external sources. Once thermalized, 
the electronic system relaxes by dissipating the energy 
in excess into the lattice and by emitting radiation. As- 
sociated to each one of the stages of the relaxation pro- 
cess, there is a characteristic time. The first one, the 
thermalization time, is an effective time determined by 
the intrinsic characteristic times of the rapid interaction 
mechanisms within the system and the second one is a 
characteristic time determined by the interaction mecha- 
nisms of the system with the surroundings. Thermaliza- 
tion and cooling processes also depend on other features 
of the excitation, for instance the excitation time, i.e., 
the time along which the carrier generation is produced 
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or the time the perturbation remains switched on. The 
thermalization and cooling times notoriously change if 
the external perturbation remains switched on over times 
longer than the e-e collision time, or if the perturbation 
switching on occurs adiabatically. 

We present in this paper a detailed study on the ultra- 
fast processes which lead to the thermalization condition 
in a photo-excited electron gas in polar semiconductors. 
We assume the semiconductor is excited by a pulsed laser. 
We define two quantities on which our discussion is based, 
namely, the thermalization time and the thermalization 
temperature We analyze how these quantities depend 
on some external variables such as the time duration of 
the laser pulse, the energy of excitation, the carrier con- 
centration, and the lattice temperature. 

II. THEORY 

To study the thermalization process of a photo- 
generated electron gas in a bulk semiconductor we use 
a set of rate equations we have developed to investi- 
gate several physical situations^, Q. These equations 
describe the evolution of the electron population in semi- 
conductors under several general conditions, and consider 
the most important interaction mechanisms, including 
screening and phonon population effects. Our theoreti- 
cal scheme is based in the following formalism. 

Let us start by defining the carrier population in a 
volume element in the (r, k)-space or /z-space, composed 
of the carrier position r and the carrier wavevector k, 

77(r,k,t) = /(r,k,i)^3 drdk (l.a) 

where /(r, k, t) is the carrier distribution function in the 
/i-space and V is the crystal volume. 

In a similar way, for the optical phonon population we 
can write 

Wi(q,f) = 9] (<l,t)j^dq , (l.b) 

where j labels the branch, q is the phonon wave vector 
and <7j (q, t) is the phonon distribution function. In the 
following we will use the index a to denote the couple 
(J, q). So that N a will denote the phonon population in 
the mode a. 

An electron in a semiconductor is characterized by its 
position r, its wavevector k, and the index of the energy 
band. An electric field introduces a preferred direction. 
By assuming an homogeneous, isotropic system, and for 
null applied electric field, the only relevant variable is 
the energy, as is the case of cubic semiconductors like 
GaAs[6l 1 1 II] - Therefore, we can write for the carrier and 
phonon populations 

/x r , . VdSde ,„ . 

^ = ^ 47r 3 |V k e(k)| < (2 - a) 

N a (t) = gj (q,t)j^dq , (2.b) 



where dS is a surface element on the surface of constant 
energy e. 

Now we wish to establish the equations which gov- 
ern the time evolution of these quantities. In order to 
do that we assume the following. Transport and op- 
tical properties of cubic semiconductors like GaAs are 
explained in terms of their band structure. The band 
structure of a cubic-model semiconductor is composed 
of one conduction band with three sets of minima, and 
three valence bands. The minima of the conduction band 
are located at the T point (k — 0), at the L points 
[k = (tt/clo, Tt/ao, 7r/ao)j being the lattice constant], 
and along the A lines k = (k, 0, 0). The tops of the va- 
lence bands are located at T point. Two of them are 
degenerate at this point and the other is separated by 
spin-spin interaction. Then, we suppose a band struc- 
ture composed of the valley T separated from the valley 
L by an energy A. This model is appropriate for the de- 
scription of photoexcited semiconductors where the tran- 
sitions involve only the center of the Brillouin zone [ill]. 
We assume also that the continuum of states of the valley 
is partitioned into a set of discrete intervals of energy Ae. 
For simplicity, and in order to have a direct reference we 
choose Ae to be the longitudinal optical (LO) phonon of 
energy Hujlo- This choice is not an essential assumption 
for the development of our procedure @. 

Instead of establishing evolution equations for the 
quantities rj(e,t), we shall set up the evolution equations 
for these quantities integrated in the interval of range 
Ae = fi^LO- So, we have 

V-yi(t) = f-f(ci,t)d- fi , (3) 

where 7 = T,L labels the valleys and i the levels of energy 
e, = «Ae; i = 0, 1,2, . . . , and 

f r> / w f VdSde 

JAe JAe 47r J |V k e(k)| 

is the number of energy states in the interval character- 
ized by set of indexes (ji) . Obviously it depends on the 
density of states D 1 (e) in the respective valley. 

We obtain the evolution equation for the carrier pop- 
ulation in the different levels in the energy space by us- 
ing the conservation of the electron number, thus we can 
write 

j~ — ^ ] (bm ~ a m) + G~fi — Ryi , (5) 
m 

where a m (b m ) is the flux out from (into) the level ji 
due to the interaction mechanism labeled by m. The 
photo-generation of carriers is accounted for G 7! ; while 
Rji denotes the recombination rate. The rate equations 
([5|) govern the evolution of the carrier system and form a 
set of non linear coupled differential equations where the 
fluxes depend on the carrier populations 
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A. Scattering mechanisms 

We now particularize our treatment to the case of po- 
lar semiconductors. However, the adaptation of the the- 
ory necessary to deal with covalent semiconductors is al- 
most direct. The evaluation of the probabilities associ- 
ated to the collision mechanisms can be done by using 
the straightforward first-order perturbation theory. The 
expressions obtained by the use of the Fermi Golden rule 
for the transition probabilities due to the different inter- 
action mechanisms, can be found somewhere else. The 
derivation of some of them and modifications of these 
expressions according with our theoretical framework, is 
almost direct, here we just discuss those details which in 
our opinion might need some clarification. Full details 
can be found in Ref. [6(. 



1. Electron-longitudinal-optical phonon interaction 
(polarization potential) 

When a carrier undergoes an electron-logitudinal- 
optical phonon interaction makes a transition to the neig- 
bour levels (i ± 1) in the same band. Transition to dif- 
ferent bands (inter valley) in which an optical phonon is 
participating is due to a different potential, the optical 
deformation potentialpj]]. For photoexcited semiconduc- 
tors we are interested in the transitions occur only at the 
center of the Brillouin zone, and these inter valley transi- 
tions are neglected [111]. The expression for the fluxes due 
to the electron-logitudinal-optical phonon interaction is 
derived in Q and are based on the matrix elements de- 
rived elsewhere, see for example II 21 . 

The fluxes can be written as[^. 1 12L ll3T] 



where 



y 7± 
op 



(0 



/m 7 e fkJho 



,7± 



±i 



N„ 



(6.a) 



(6.b) 



1 1 
- ± - 

2 2 



1 fruJLO 



±lTV / lT ^LO 



(7) 



Sod and £ s are the static and optical dielectric constants, 
respectively. m 7 is the effective mass in the 7 valley, 
e is the electron charge and h is the Planck constant. 
The upper (lower) sign is for emission (absorption). The 
screening effects in the e-ph interaction are included in 
the factor S^O' which in the random- phase approxima- 
tion is given by[13j 
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where TV and iV 7 are the total and 7 valley carrier con- 
centration respectively, is the threshold value for the 
concentration in the 7 valley at which the screening be- 
comes important. This critical value for the carrier con- 
centration is given by 



N c = 

7 



m 7 (huj LO ) 3 
3 3 / 2 8ire 2 h 2 k B T c 



where T D is the electronic temperature and fee is the 
Boltzmann's constant. 



2. Electron-electron interaction 

The electron-electron interaction gives the nonlinear 
character of the Eqs. (5) and is one of the most difficult 
interactions to take into account. We adopt the Debye- 
Hiickel screened potential to describe the e-e interaction. 
We see that the scattering processes in which the mag- 
nitude of exchanged momentum is small are the most 
likely, because de probability diminishes as q~ A . Also, 
on the average, an electron in the valley T exchange ap- 
proximately fuvLO of energy Q. Therefore, in our energy 
levels scheme, carriers make transitions to the neighbour 
levels i ± 1 in the same band due to the e-e scattering. 
For more details see Refs. Jp. H^|. 

The fluxes are given by[6| 



1 - 



d-fi-i 



1 - 



d. 



d-yi+l 
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1 - 



I7V+1 



d-y'i' - 



(8.a) 



Z~ 



■y'i' 



1 - 



d-y'i'+l 



'y'i' 



1 



"■y'i' 



(8.b) 



where Pauli exclusion principle has been taken into ac- 
count and 



Z'< = 



e 2 ' irm^k-QTc 



1 



1 



N 
NZ 



(9) 



The square brackets factor takes into account the screen- 
ing effects. These become important when the carrier 
concentration N reaches a critical value N^ c given by 



N 1 = 



4?ti 7 £ 00 (fceTe) 2 



The expression for the total probability, Eq. ©, is 
an heuristic useful expression, that allows us to deter- 
mine in an easy way the ranges of carrier concentration 
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and electronic temperature, in which the energy exchange 
through e-e scattering is the dominant mechanism in the 
kinetics of the system. 



B. Generation and recombination 

The dynamics of carriers under intense laser irradia- 
tion has been described by Ferry [15j. Here we adopt the 
following procedure to describe the generation and re- 
combination processes. The rate equation which governs 
the effect of the generation and recombination processes 
on the CDF can be written as 



dt 



= G 



n 



R 



n 



(10) 



the first r.h.s term represents the photo-generation of 
carriers and the second one represents the recombina- 
tion processes. In direct gap semiconductors the photo 
excitation involve the top of the valence band and the 
bottom of the conduction band which are at the T point. 
For indirect gap semiconductors it is necesary the par- 
ticipation of a phonon in order to conserve momentum. 
In this report we are interested in direct gap semicon- 
ductors, therefore, we shall assume only generation of 
carriers to the T valley, therefore, the generation term 
can be written aslal 



where i p Ae is the excitation energy (measured from the 
bottom of the conduction band) . G p (t) is proportional to 
the rate of generation, i.e., the number of excited carriers 
per unit time given by 



G p (t) = G*G(t) 



AdhujL 



G(t) 



where Pl is the laser power, Swl is the photon energy, d is 
the penetration length, and A is the area of illumination. 
The dimensionless function G(t) is conveniently chosen 
in order to model a given experimental condition. For 
instance, for the case of a steady state G{t) = 1. The 
factor a is given by /// where / is the hole distribution, 
and in general is a function which depends upon e and 
t. In steady state processes a accounts for the degree of 
compensation in the material. 

The rate of recombination can be expressed as R = 
i? R + i? N , i.e., it is the addition of the radiative and non 
radiative components of the recombination. 

In this work we assume Rn — wrjn where w is defined 

as 



1 



W 



(12) 



and T rec has to be calculated for the pertinent kind of 
recombination. 



C. Phonon population effects 

We have mentioned that a high population of phonons 
can produce some important effects and these effects are 
more notorious in the case of the LO phonons [HI . How- 
ever, the number of LO modes is limited by the mag- 
nitude of the electron wave vector. These limits can be 
easily obtained by applying the energy and momentum 
conservation conditions to a transition in which an elec- 
tron with an energy near the maximum energy e max ab- 
sorbs a LO phonon. In this way one obtains 



(7min 



7 c max 



and 



7 wnax 



l 



l 



Of course in this energy scheme it is not possible to know 
the wave vector of the absorbed phonon. In order to 
describe phonon population effects we need to link the 
carrier system to the phonon population. To this end we 
use the following evolution equation for the whole phonon 
population TVlo 



dN LO 
dt 



1 



Nil 



TLO 



(13) 



where U is the number of permitted modes per unit vol- 
ume, 



U = 



<?max 



dq 
(2VF 



./V£q is the equilibrium total LO phonon population and 
tlo is its time life. 

The set of equations §5§ and (fT3"]) with the respective 
expressions for the fluxes are the basis of this model. 
These equations can be easily extended to represent a 
more general situation. For instance, the hypothesis of 
parabolic valleys can be changed to allow the descrip- 
tion with a more realistic band structure. The solution 
to these rate equations can be considerably simplified 
by the use of effective collision frequencies (ECF). These 
frequencies are the average over the band of the scat- 
tering frequencies defined for each of the energy levels 
and which are the factors appearing in the fluxes. We 
would like to notice here that the information about the 
band structure and other symmetries of the system is 
contained in the ECFs. Thus, this approach can be used 
to study ultra-fast phenomena in quantum wells, super- 
lattices and heterostructures, as well as bulk semiconduc- 
tors systems @, [13] • 
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III. NORMALIZED RATE EQUATIONS 

For not very high carrier concentration we can neglect 
the inter valley transitions. The threshold in which we 
can neglect inter valley transitions depends on the ma- 
terial, for GaAs this threshold is 1 x 10 18 cm -3 [Hi. If 
additionally one assumes a non-degenerate electron gas 
and use the ECFs as defined before, the rate equations 
([5]) and (fT3|) become notoriously simplified. The use of 
the ECFs in our description is in fact justified by the 
smooth behavior of the scattering probabilities as func- 
tion of the energy. Hence, under these conditions the rate 
equations can be cast into 

+ ZN max x{Xi+i ~ 2 Xt + Xi-i) 

+ ZN max xo(Xi ~Xi-i) 

+ 9 P Ki v -wXi , ( 14 -a) 

for i ^ and for i — we have 

dxo + 

+ ZN max x(Xi ~ Xo) + ZN^^xoXo 

+ 9 P So,i p - wxo • (14.b) 

For the phonon population in excess, Nf^ — Nlo — N^q, 
we have 

JAfCXC 1 

= -K(x-Xo)-^x]-tN^ (14.c) 

where the populations have been normalized to the max- 
imum reachable carrier concentration -/V max , i.e., Xi = 
W^max; X = SjXi; u = U/N max ; and g p = G p /N ma *. 
The ECF u£ has just the dominant term of the e-LO 
phonon interaction, and £ = 1/tlo Notice the differ- 
ences between the rate equation for i = and i 7^ 0. 



IV. RESULTS AND DISCUSSION 

For the sake of brevity, from here on, we will refer as 
carrier distribution function to the set of values {xi} of 
the carrier population at the different energy intervals on 
the conduction band. The link of this distribution with 
the actual out of equilibrium CDF is given by the expres- 
sion ([3]). We start our discussion by defining two physical 
quantities inherent to the time evolution of the photo- 
generated carrier population. More specifically, we wish 
to characterize by means of these physical quantities, the 
stage in which the system reaches the thermalization con- 
dition. These quantities are the thermalization time t* 
and the thermalization temperature T*pJ. Our defini- 
tion of these quantities intend to be phenomenologically 
amenable. Thus, t* is defined as the time interval, mea- 
sured from the begining of the laser pulse that the carrier 



system requires to reach a distribution shape which can 
be fitted by means of a single exponential function, 

Xl = AexpH'Ae/fc B T e *) (15) 

Here A is a normalization constant. The quantity T*, 
which makes the fitting possible, defines the thermaliza- 
tion temperature. This last definition closely resembles 
the way in which the carrier temperature of a hot elec- 
tron system is experimentally determined We proceed 
by numerically integrating Eqs. (114. al e) under different 
physical conditions and seeking the effects of the external 
variables on t* and T*. Starting at the pulse begining, 
the carrier population generated in the conduction band 
evolves due to the collision mechanisms according to the 
rate equations (|14.al -c). At each step of the numerical 
integration we are able to calculate the CDF and the 
corresponding values of the ECFs. In the initial steps 
of the numerical integration, in general the CDF differs 
clearly of the shape of a thermalized distribution, i.e., 
a Boltzmann factor. Within subsequent iterations the 
CDF gradually addopts a decreasing shape, which even- 
tually, at t = t*, admits a fitting by means of a simple 
decreasing exponential function. The exponent in that 
event is inversely proportional to T* [J H, @] • This is 
the way in which we proceed to evaluate t* and T*. 

In order to make concrete calculations we consider 
the well known values of the electronic band structure, 
phonon dispersion relations and material parameters of 
GaAs. For example, fiw^o = 36 meV and the energy ex- 
tent of the conduction band at the T point is 1 eV. This 
means that the necesary number of energy levels is 28. 
Other GaAs material parameters can be found in [|| [ll| 
and references there in. At t — in an empty conduction 
band a laser pulse of duration t p injects g p electrons per 
unit time with an energy i p in units of Ae and above of 
the bottom of the band. 

Before starting our analysis of the thermalization pro- 
cess we discuss briefly the physical suitability of our car- 
rier temperature definition. By means of expression (|15p 
we have defined at time t* the thermalization temper- 
ature T* . At subsequent times one also might use this 
procedure to calculate the carrier temperature T e (t). We 
have studied the time evolution of this quantity under 
various different physical conditions. Our results show 
a good agreement with experimental data@. Firstly we 
analyze the effect of the pulse duration on the time evolu- 
tion of the carrier temperature T e , as well as on the main 
interaction frequencies, i.e., e-e, e-ph emission, and e-ph 
absorption; z, v + , and v~ respectively. In Fig. 1(a) we 
show T e (t) for an instantaneous pulse t p = ps which ex- 
cites = 10 17 cm -3 electrons at an energy level i p = 4. 
We have also plotted the behavior of the kinetic tempera- 
ture defined by 2(e)/3fcs. We show the behavior of these 
quantities since the pulse starts, until times well above 
the thermalization condition is reached. We assume in 
this calculation a lattice temperature Tl = 10 K and 
take from the literature a commonly used value of the 
damping constant for the LO-phonons r^o = 12 ps. No- 
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FIG. 1: Time evolution of the effective carrier temperature T e 
and the kinetic temperature obtained from the mean energy 
in units of K, 2{e)/3fcs. The time evolution of the main ECFs 
are shown in the lower part. The pulse duration is (a) t p — 
ps and (b) t p = 1 ps. 



tice that the resulting thcrmalization time is t* — 0.26 ps 
(indicated by an arrow in the figure). In the lower part 
of Fig. 1 (a) appears the time evolution of the main ECFs 
during this relaxation process. In Fig. 1(b) we show these 
quantities for a pulse of finite duration t p = 1 ps. In this 
case we obtain t* = 0.65 ps (also indicated by an arrow in 
the figure) . From the comparison of these figures we may 
conclude that the pulse duration has a clear influence on 
the general characteristics of these quantities, namely, 
the pulse duration modifies t* and T*. Thus, this is one 
of the external variables which determines the interaction 
mechanisms and their relative importance in the thermal- 
ization process of the CDF. The evolution of the main 
ECFs provide detailed information, at every stage of the 
relaxation process, about the relative incidence of the 
main collision mechanisms in the kinetics of the carrier 
system. We wish to recall here that we are considering in 
our treatment screening and phonon population effects. 
We notice that the e-e ECF varies rather slightly along 
the period of time considered, but the ECFs correspond- 
ing to e-ph interaction exhibit more pronounced varia- 
tions. Notice the semi-logarithmic scale of the graphics. 
In Fig. 1(b) one can observe a clear abrupt change in the 
time evolution of the electronic temperature at the time 
when the pulse ends. We now restrict our analysis to the 
thcrmalization process and, on the basis of the behavior 
of the ECFs in this stage, we discuss the effect of t p and 
other external parameters on T* and t*. 




FIG. 2: Temperature and time of thermalization vs pulse du- 
ration for two values of the excitation energy level. Contin- 
uous lines correspond to i p = 4 and dashed lines to i p — 10. 
The carrier concentration is N = 10 17 cm -3 . The scale for t* 
is at the right axis of the graph. 



Fig. 2 shows the changes which the pulse duration in- 
duces on T* and t* . We consider here a carrier population 
N = 10 17 cm~ 3 , a lattice temperature Tl = 10 K, and 
two energy levels of excitation i p = 4 and i p = 10. Notice 
that the scale for t* appears on the right hand side of the 
graphs. We observe that both T* and t* reach higher 
values for i p = 10 than they do for i p — 4, this is so all 
along the interval of t p . This result is easily explained if 
we realize that for i p = 10 the carrier system receives a 
higher excess energy than it does for i p = 4, it leads the 
carrier system to reach, comparatively, a higher value of 
T* and due to the also relatively larger number of en- 
ergy states accesible to the carriers, the thermalization 
condition requires a longer time. The other aspect of the 
figures worthy of mention is the difference between the 
behavior of T* and t* in the region of low values of t p , 
for both i p — 4 and i p = 10. While t* increases with 
t p up to reach a maximum value and then decreases, T* 
decreases monotonously. These last asymptotic behav- 
ior would correspond to the values of time and temper- 
ature of thermalization that the system would reach in 
the event of a CW laser mode experiment in which the 
same rate of generation is kept during a long time, i.e., 
t p — > oo. Again here the larger energy in excess received 
by the carriers and the large energy of states of the sys- 
tem with i p = 10 explain the corresponding larger value 
of T* and t* , in comparison to the respective results for 
i p = 4. 

The carrier population at the lowest level in the con- 
duction band governed by Eq. (14. c) has an important 
role in the thermalization process [l4j. The shoulder 
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Lattice Temperature [K] 



FIG. 3: T e * and t* as a function of the lattice temperature 
Tl for N = 10 17 cm -3 and two pulse durations t p — 1 ps 
(continuous lines) and t p = 10 ps (dashed lines). 



FIG. 4: T e * and t* vs energy level of excitation i p for the same 
pulses as Fig. 3. 



shown in the behavior of t* as a function of the pulse 
duration t p , is caused predominantly by two aspects of 
the carrier kinetics. One of them is the time necesary 
to form in the lowest level a carrier population, which 
in the thermalized condition must exponentially decay 
for increasing energies. The other aspect is the rate of 
carrier generation at the level of energy i p , which in the 
thermalized condition should be small enough in order 
that this generation at each step becomes included in 
the error tolerance in the fitting. 

Fig. 3 shows how the values of t* and T* change with 
the lattice temperature Tl- We consider here two dif- 
ferent pulse durations t p = 1 ps (continuous line) and 
t p = 10 ps (dashed line), an excitation energy i p = 4 and 
a carrier concentration TV = 10 17 cm -3 . We observe that 
the lattice temperature has a rather mild effect on the 
thermalization process. According to the behavior of the 
ECFs during the thermalization process (Fig. 1), the ma- 
jor influence must come from phonon absorption events, 
whose contribution is larger for higher Tl values, this is 
so, because it propitiates a high LO-phonon population 
in such a way that, comparatively to the phonon emis- 
sion, phonon absorption increases its participation in the 
kinetics of the carrier system. 

In Fig. 4, for the same pulses as in Fig. 3, we have de- 
picted our results for t* and T* as a function of the energy 
of excitation i p . The dashed curves correspond to t p = 10 
ps. The non-linear nature of the e-e interaction is rev- 
eled in these results. We observe that a general behavior 
of both quantities is that, they increase with increasing 
values of the excitation energy, however, this behavior is 
non-monotonic. The dependence of the carrier tempera- 
ture upon the carrier concentration and upon the energy 



of excitation have been studied since the 70's decade. 
We have studied this dependences theoretically. In par- 
ticular we have found that strong changes in the CDF, 
associated with the non-linear nature of the e-e interac- 
tion, can be induced by varying the excitation energy [l8j. 
The abrupt change in T* and t* about i p — 14 is a man- 
ifestation of this sensitive dependence. There has been 
in the literature some discussion regarding to this point. 
In particular, this phenomenon has been analized as a 
phase transition like behavior of the carrier population 
in the lowest level of energy 14 1 . 

In Fig. 5 we show the dependence of T* and t* on the 
carrier concentration. We also consider here two pulses 
of duration t p = 1 ps, continuous lines, and t p = 10 
ps, dashed lines, an energy level of excitation i p = 4, 
and a lattice temperature Tl — 10 K. The non- linear 
dependence of t* and T* upon the carrier concentration 
is clearly exhibited. Notice that for carrier concentra- 
tions larger than 10 17 cm -3 the screening effects begin 
to be noticeable. Screening turns the e-e interaction less 
effective in thermalizing the CDF and reduces the rate 
at which e-ph scattering takes out energy from the elec- 
tronic system, in this way, both T* and t* increases with 
the carrier concentration. 



V. COMMENTS AND REMARKS 

The criterion we have applied here to define T* and 
t* by means of the expression (|15[) may require some 
improvement. We think that the least squares fitting 
to an exponential function, although it resembles an ex- 
perimental procedure, the included inherent error could 
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FIG. 5: Dependence of T e * and t* on the carrier concentration. 
Energy of excitation i p — 4 and pulses of t p = 1 ps (continuous 
lines) and t p = 10 ps (dashed lines). 



be the origin of some of the uneven behavior we ob- 
serve in various of our results. However, the agreement 
of our results in the description of the cooling and the 
steady state processes in hot electron systems @, with 
the experimental data, provides support to our theoret- 
ical model to describe the ultra-fast phenomena which 
lead a photo-generated carrier population to reach the 
thermalized condition. The reported experimental and 
theoretical results on the subject of the kinetics of the 
thermalization^ also support our definition. One of the 



advantages of our "kinetic" approach is that, due to the 
fact that ECFs depend only on the band structure of the 
system and consequently on the system dimensionality, 
it can be applied to study some ultra-fast phenomena in 
systems of reduced dimensionality and also systems with 
a small number of particles. In fact we have applied these 
theoretical framework to study some transport phenom- 
ena in mesoscopic semiconductor heterostructures [1 7j . 
An additional interesting characteristic of this theoreti- 
cal procedure is that, the required numerical calculations 
are not at all expensive. All the results we present in this 
paper can be obtained in a few minutes in an ordinary 
PC. In our analysis we have focused our attention to the 
role that the main interaction mechanisms play, however, 
the inclusion of some other scattering processes in the ki- 
netics of the system is a simple task in this theoretical 
scheme @. 

In conclusion, we have presented a detailed analysis 
of the thermalization process in terms of the relevant 
external parameters, i.e., laser pulse duration, energy of 
the photoexcitation, intensity of the photoexcitation, and 
lattice temperature. We defined two physical parameters 
inherent to the time evolution of the system on which the 
analysis of thermalization have been done, the thermal- 
ization time and temperature. We have found that the 
lattice temperature has a neglegible influence on T* and 
t* . The other three external parameters clearly influence 
the thermalization process. 
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